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Abstract 



We study the contribution of the squark flavor mixing from the LR(RL) component of 
the squark mass matrices to the direct CP violation of the b — > S7 decay and the CP- 
violating asymmetry in the non-leptonic decays of B mesons. The magnitude of the 
LR(RL) component is constrained by the branching ratio and the direct CP violation of 
b — > sj. We predict the correlation of the CP asymmetries among A b ^, sl , S^Ks an d S v i K o 
of the B decays. The precise data of these CP violations will give us the crucial test for 
our framework of the squark flavor mixing. 
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1 Introduction 



New physics are expected to be observed at the LHC experiments. Although new particles 
have not yet discovered, LHCb has reported new data of the CP violation of B mesons and 
the branching ratios of rare B decays. New physics are also expected to be found in the B 
meson decays. 

The CP violation in the K and B^ mesons has been successfully understood within the 
framework of the standard model (SM), so called Kobayashi-Maskawa (KM) model [1]. The 
source of the CP violation is the KM phase in the quark sector with three families. However, 
there could be new sources of the CP violation if the SM is extended to the supersymmetric 
(SUSY) models. The CP-violating phases appear in soft scalar mass matrices. These phases 
contribute to flavor changing neutral currents with the CP violation. Therefore, we expect 
the SUSY contribution in the CP-violating phenomena in the B meson decays. 

The typical contribution of SUSY is the gluino-squark mediated flavor changing process 
[2]-[TT]. In our previous paper [12], we have already discussed the effect of the squark flavor 
mixing on the CP violation in the non-leptonic decays of B% and B® taking account of the 
recent LHCb experimental data. We have found the deviation from the SM predictions in 
the asymmetries of the penguin dominated decays B® — > <pKs and B® — > rj'K . In that 
framework of the SUSY contribution, we assume that LR and RL components of the squark 
mass matrices are neglected. The LL and RR components of squark mass matrices contribute 
considerably to the penguin processes for the case of large \i tan j3, O(10 TeV). However, if LR 
and RL components of squark mass matrices dominate the penguin decays, the asymmetries 
of B® — > (j)Ks and B® — > f]'K° are deviated from the SM predictions even for the case of 
the smaller /^tan/3, (9(1 TeV). Then, these contributions of the new physics are correlated 
with the direct CP violation of the b — > 57 decay. In this paper, we present the numerical 
analyses in the case that LR and RL components of squark mass matrices dominate the 
penguin decays. In this case, the LR(RL) components do not contribute to the dispersive 
part Mf 2 of B q — B q (q = d, s) mixing. 

In section 2, we summarize the effect of new physics in the CP violations of the neutral B 
mesons including the recent experimental data. In section 3, we discuss the our framework 
of the squark flavor mixing in the CP violation of B mesons. We also discuss the constraints 
from the direct CP violation in the b — > 57 process. In section 4, we show the numerical 
result of the CP violation in the B mesons. Section 5 is devoted to the summary. 



2 New physics of CP violation in B mesons 

Let us discuss the effect of new physics in the non-leptonic decays of B mesons. The contri- 
bution of new physics to the dispersive part Mf 2 (q = d, s) is parameterized as 

M\ 2 = Mff u + Mli USY = Mlf M (l + V 2 " 9 ) > (Q = d > s) (1) 

where M^ SUSY is the SUSY contribution, and the SM contribution M^ SM is given as [T3~| 

Mli m = ^^MUV tb V:^ B S (x t )f B B q . (2) 
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The time dependent CP asymmetry decaying into the final state / is defined as [H] 



_ 2ImA ; 



where 



x _q. q_ I Mil - |r; 2 _ _ A(B° q -> f) 

In the .8° — >• J/ipKs decay, Xj/^k s is given in terms of the new physics parameters hd and 
ad as 

Aj/^ s = ~e-^ d , <p d = 2(3 d + arg(l + h d e 2 ^), (5) 

by putting \p\ = 1 and q/p ~ a/ Mf;/Mf 2 , where the phase ^ is given in the SM. The 
CKMfitter provided the allowed region of hd and ad, where the central values are [T5l IT6] 



hd ^ 0.3, crrf ~ 1.8 rad. (6) 

Since penguin processes are dominant in the case of / = (j)Ks,T]'K , the loop induced new 
physics could contribute considerably on the CP violation of those decays. Then, Sf is not 
any more same as Sj/^k s if new physics leads to \p\ ^ 1. Those predictions provide us good 
tests for new physics. 

In the B° s — > J/ip<f) decay, we parametrize as 

A W = e"^, 4>s = -2/3, + arg(l + h s e 2i °*), (7) 

where /3 S is given in the SM. Recently the LHCb presented the observed CP-violating phase 
(f> s in — > J /ipit + it~ decay using about 1 fb" 1 of data [17] . This result leads to 

<j> s = -0.019™+^ rad, (8) 
which is consistent with the SM prediction [T5] 

tfmsM = _ 2/3g = _ .0363 ± 0.0017 rad. (9) 

Taking account of these data, the CKMfitter has presented the allowed values of h s and a s 
[T5| [T6j . The allowed region is rather large including zero values. In order to investigate 
possible contribution of new physics, we take the central values 

h s = 0.1, a s = 0.9 - 2.2 rad, (10) 

as a typical parameter set in our work. 

We remark on numerical inputs of phases <f) d an d <f) s in our calculation. The phase <pd is 
derived from the observed value Sf = 0.671 ± 0.023 in B% — > J /ipKs PS] as seen Eqs.flB]) 
and (EJ). On the other hand, we use the SM value of f3 s and the values of the new physics 
parameters, h s and a s in Eq.lflOl to estimate (p s = —2(3 S + arg(l + h s e 2lC7s ). We do not use 
the observed value of <p s in B° s — > J/i/j4> because of the large experimental error in Eq.([8]). 
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Since the —> J/ifjKs process occurs at the tree level in SM, the CP- violating asymmetry 
originates from Mf 2 . Although the B% — > <pKs and B% — > rj'K decays are penguin dominant 
ones, their asymmetries also come from Mf 2 . Then, asymmetries of B® — > J/tpKs, B® — > <pKs 
and B® — > r]'K are expected to be same magnitude in SM. 

On the other hand, if the squark flavor mixing contributes to the decay at the one- 
loop level, its magnitude could be comparable to the SM penguin one in B® — > (pKs and 
Bd rfK , but it is tiny in B® — > J/ipKs- Endo, Mishima and Yamaguchi proposed the 
possibility to find the SUSY contribution in these asymmetries [19]. The present data suggest 
the deviation from SM in these time dependent asymmetries of B^ decays such as, 

Sj/ipk s = 0.671 ± 0.023, S^ Ks = 0.39 ± 0.17, S v , K o = 0.60 ± 0.07, (11) 

however, precise data are required to justify the new physics contribution. 

New physics contribute to the b — > 57 process. The observed b — > 57 branching ratio (BR) 
is (3.60±0.23) x 10" 4 [18], on the other hand the SM prediction is given as (3.15±0.23) x 10~ 4 
at 0(a 2 s ) [201 [21]. Therefore, the contribution of new physics should be suppressed. New 
physics are also constrained by the direct CP violation 

A b*r = r(£ -» X sl ) - T(B -» Xry) 

CP r(B^x sl ) + r(B^x sl y { ' 

Since SM prediction A^ S7 — 0.005 is tiny [22], new physics may appear in this CP asym- 
metry. The present data A b ^, sl = —0.008 ± 0.029 [18] has rather large error bar, and so the 
constraint of new physics is not so severe. However, improved data will provide the crucial 
test for new physics. 



3 Squark flavor mixing and CP violations of B mesons 

Let us consider the flavor structure of squarks in order to estimate the CP- violating asymme- 
tries of B meson decays. We take the most popular anzatz, which is to postulate a degenerate 
SUSY breaking mass spectrum for down-type squarks. Then, in the super-CKM basis, we can 
parametrize the soft scalar masses squared M~ , M ~ , M~ , and M~ for the down-type 
squarks as follows: 

(#% (#% (#% 



+ (sr)n OCT 

M2 d=<\ 1+or 



23 



d RR 1 



d 

LR\ 
d JH 



X«CTsX 




122 



M L = ( M !J = m l\^ R ^ 0CT*>l> (13) 
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where rriq is the average squark mass, and (8 d L )ij, (^d R )ij, (8d L )iji an d (^d R )ij are called as 
the mass insertion (MI) parameters. The MI parameters are supposed to be much smaller 
than 1. 

The contribution of the gluino-squark box diagram to the dispersive part of the effective 
Hamiltonian for the B q -B q mixing is written as 



m1 ,susy = A , ^ {{5 l L% + {6SR)3 j + m ^ U5 R R)i 



+ At {(5 L d % + {5*%} + Al{5 L d R U5 d 



RL\ 



ii.l 



(14) 



where 
A\ 



A\ 



A\ 



216m? 3 



MbJL, 



384 



-132 



A 2 = 24xf 6 (x) + QQf 6 (x), 




(15) 



Here, we use x = rn^/m 2 -, where m~ g is the gluino mass. For the cases of q = d and q = s, we 
take = (1,3) and = (2,3), respectively, where m\ = m d , m 2 = m s and m 3 = m b . 
The loop functions f%{x) and f%{x) are shown in 



For the case of x ~ 1, we get A 2 



■1, A\ ~ 30, A| ~ 



■10 and A; ~ 10. Therefore, each 



term at the r.h.s. of Eq.( fl4|) may contribute to Mf 2 SUSY comparably. However, magnitudes 
of (S d R )ij and (8 d L )ij are constrained severely by the b — > 57 decay as discussed later. 

The squark flavor mixing can be tested in the CP-violating asymmetries in the neutral B 
meson decays. Let us present the framework of these calculations. The effective Hamiltonian 
for AB = 1 process is defined as 



H. 



eff 



4G F 



_q'=u,c i=l,2 

where the local operators are given as 
>(<?') 



E E c *°^ - W E + Si *) 



i=3-6,77,8G 



(16) 



0} 9J = (s a7 M^)(?>W a ), 
3 = (5 a7M PL&a) E(^7^L^) 

1 

5 = {Sal^Lba) E(^7 M -PR^; 



or 



(SaJpPLq^tfp-fPLbp). 

4 = (saipPLbp) ^{q^P L qc 



o 



167r : 



-,m h s ot a tlv P R b ol F IJtu , O 



6 = i.S a ^^P L bp) E(^7 M ^a) 



7 8G 



167T 



; m b s a a^P R T°bpGl 



(17) 



where Pr = (1 + 7s)/2, P^ = (1 — 7s)/2, and a and /3 are color indices, and g is taken to be 
it, d, s, c. Here, Cj's Cj's are the Wilson coefficients, and CVs are the operators by replacing 
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L(R) with R(L) in Oj. In this paper, C, includes both SM contribution and gluino one, such 
as C{ = Cf M + Cf, where Cf M is given in Ref. [25] and Cf is presented as follows [26] : 



CI 



LL\ 



c 9 
cf 



"WW 1 
v/2V 2 



ts" u q 



~ 4G F V tb V t * s mf d m 



LL\ 



V2ai 



;LL\ 



~ AG F V tb V t * s mj 



2^d m 



- l -B x {x) - ~B 2 (x) - iPi(x) - l -P 2 (x) 

- 7 -B x {x) + l -B 2 {x) + ipx(x) + |p 2 (x) 

10 1 1 1 

—Bxtx) + -S a (x) - —Pilx) - -P 2 (x) 

-\b x {x) + ? -B 2 (x) + ipx(x) + |p 2 (x) 



C; 



\/2< 



7o 



6C F ^m? 



m| 3 



(<^) 23 I ~M 3 (z) -^tan/3^M a (x) ) + (^) 23 ^jM 1 (x 



mj 3 



n 



>/2. 



<s<v 



2G F T4 6 y t >? 



(^ L ) 23 <; ( -m 3 (x) + 3m 4 (x; 



,m 5 / 1 



//tan/3— § -M a (x) + 3M 6 (i; 



(0 23 ^(k(x) + 3M 2 (x; 
m 6 V3 



Here the double mass insertion is included in C| and Cf G . The Wilson coefficients Cf s 
are obtained by replacing L(R) with R(L) in Cf s. The loop functions, which we use in our 
calculations, are presented in our previous paper [T2] , 

The CP- violating asymmetries Sf in Eq. (J3J are calculated by using A/, which is given 
for -> and 5° ->■ T)'K° as follows: 



j= 3 -6,77,8G 



i= 3 -6,77,8G 



(19) 



where (O t > is the abbreviation of </|Oi|^J>. It is noticed that (^A'slO^) = (0^10^5°) 

and (rj'K°\Oi\Bj) = —(r]'K°\Oi\B^) because of the parity of the final state. We have also A/ 
for B® — > 66 as follow: 



-i*. ^=3-6,7 7 ,8G 

i= 3 -6,77,8G 



(20) 



with (#10,1^°) = -(#10*15°). 
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In these non-leptonic decays, the Cg G (Oso) dominates these amplitude, but small contri- 
butions from other Wilson coefficients are also taken account in our calculations. Therefore, 
we estimate each hadronic matrix elements by using the factorization relations in Ref. |27j . 

Let us discuss the each contribution of the mass insertion parameters to C| G in Eq. ffTgj) . 
Taking account that the loop functions Mi(x) are of same order and m g ~ rriq, the ratio of 
LL component and LR one is (6^)23 x //tan (3 /rriq to (6^)23 x mq/m b . If 0(/itan/3) ~ 
0(rriq) and rriq > 1 TeV, the LR component may contribute significantly to C| G due to 
the enhancement factor mq/m b = O(10 2 ). For example, in the case of (6^)23 = 10~ 2 and 



(5^)23 = 10 3 , the LR component dominate Cf G , while it is minor in M\ 2 



23 

SUSY 



as seen m 



Eq. (jM|) . This situation is also kept in the b — > 57 decay. 

The b — > S7 decay is a typical one to investigate new physics, 
given as 

BR(B X sl ) 



\v; s v tb 


2 Qa 






2 





where 



BR(B X c ev e 



f(z) = 1 - 8z + 8z 3 - z A - I2z 2 lnz 



.^fcff|2 
°7-7 I J 



The branching ratio is 



(21) 



mi 



(22) 



Here, includes both contributions from the SM and the gluino-squark flavor mixing G 



As seen in Eq. ffT8j) . both and C| G have the similar dependence of (6^ )23. Therefore, we 
should discuss carefully the contribution from (5^)23 in our numerical calculations. 

We can discuss the direct CP violation A 1 ^ 111 in the b — > 57 decay, which is given as [22] 

T(B -» X sl ) - T(B -» X s - 7 ) 
T(B -)■ X s7 ) + r(B -> X s - 7 ) 

r4D 87 
^Im[C 2 C; 7 ]-^Kz) + ^,5)]lm 



^CP 



£ 7 >(l-5)£2 



IC- 



■7 7 1 

- ^Im[C 8G C^] + *)Im 



V us V ub 

v t * s v tb 



GoG 



2^77 



u-6 



CoCi 



where v(z) and 6(^,5) are explicity given in [22], and C 4 eff includes both the SM and SUSY 
contributions. Although the experimental data has still large error bar, we can discuss the 
SUSY contribution to A^ 111 . 

Let us set up the framework of our calculations. Suppose that /xtan/3 is at most 0(l)TeV. 
Then, magnitudes of (6^)23 and (64)23 are constrained by M( 2 as seen in Eq. (|T4|) . Taking 
account of h s = 0.1 in Eq. flTU]) . we obtain | (5^)23! — \(6d R )2s\ — 0.02 in our previous work 
[12]. Then, these contributions to C|L and C| G are minor. 

On the other hand, (5^)23 and (6^)23 are severely constrained by and C§ G . We 
show the constraint for (64 ) 
we suppose 1(5^) 23 1 = |(^ L 



23 and (6^)23 in our following calculations. In our convenience, 
) 23 1 - Then, we can parametrize these parameters as follows: 



,2ir9 



LR 



(^)23 = |(^)23|e^, 

where 8^ and 9^ are taken in the region [0 
analyses in the next section. 



(O23 



RL 



|(023|e 2 ^ 



(23) 



71 



By using this set up, we show numerical 
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4 Numerical analyses 



In this section, we show the numerical analyses of the CP- violation in the B mesons. In our 
following numerical calculations, we fix the squark mass and the gluino mass as 

rrig = 1000 GeV, m~ g = 1500 GeV, (24) 

which are consistent with recent lower bound of these masses at LHC [28]. We use relevant 
parameters as given in [12] to estimate the SM contribution. 

At first, we discuss the b — > 57 decay. The observed b — > 57 branching ratio is (3.60 ± 
0.23) x 10 -4 p2], on the other hand the SM prediction is given as (3.15 ± 0.23) x 10~ 4 at 
0(a 2 s ) [20], ETJ . Since /xtan/3 is supposed to be lower than 0(1 TeV), the contribution of 
(5^)23 — (5^)23 is negligibly small. The contribution of (5^)23 becomes important through 
the interference with the SM component in the decay amplitude. On the other hand, since 
(^ i )23 does not interfere with the SM component, its contribution is minor. The branching 
ratio gives the constraint for the magnitude of (5^)23- The direct CP violation of the b — > S7 
is also useful to constraint (S^ R ) 2 3- We show the A 1 ^ 1 versus 1(5^)231 in Figure 1, where 
the upper and lower bounds of the experimental data with 90% C.L. are denoted red lines, 
and the predicted value of the SM is shown by the green line as the eye guide. As far as 
1(^^)23! < 10~ 3 , the predicted value is within the experimental allowed region. 

In Figure 2, we show the 1(^^)23! dependence of the branching ratio taking accont the 
constraint of A 1 ^ 1 as seen in Figure 1. Here, the allowed region at | ((5^^) 23 1 = is the SM 
prediction. As the magnitude of (5^)23 increases, the predicted region of the branching ratio 
splits into the larger region and smaller one. The excluded region between two regions is due 
to the constraint of A b ^, SJ . Then, the predicted branching ratio becomes inconsistent with 
the experimental data at 1(^^)23! > 5.5 x 10~ 3 . 

In order to see the role of the phase 9^, we show A b ^ sl versus 6^ for | (5^)23! = 
10~ 3 (blue) and |(<5^ R )23| = 10 _4 (orange) in Figure 3. The pink horizontal lines denote the 




l(^ K )23l 



Figure 1: The direct CP violation A^ 1 
versus 1(5^)231, where the green line de- 
notes the SM prediction and red lines de- 
note the upper and lower bounds of the 
experimental data with 90% C.L.. 



_ 4.oxio- 4 
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g 3.0X10- 4 

2.0X10" 4 





- iur''t'"' : * ■'. Vr*' ; 













2.0xl0" 3 4.0xl0" 3 6.0xl0" 3 



l(^ R )23l 

Figure 2: The predicted branching ratio of 
b — > 57 versus 1(5^)23], where the exper- 
imental constraint of A b ( ^, s ' y is taken ac- 
count. Predicted value at \(8% R )23\ = is 
the SM one. 
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experimental upper and lower bounds at la level. As seen in this figure, we find that the 
reduction of the experimental error-bar will constrain the SUSY phase 9^ severely. 

In Figure 4, we plot the allowed region on the | (5^)23! — plane by putting the 
experimental data at 90% C.L. of the branching ratio and the direct CP violation . 
The allowed region of | ( < 5 ( f' jfe ) 23 1 is cut at 5.5 x 10~ 3 , where 9%^ is tuned around ir/2. Around 
7r/4 and 3"7r/4, ^4^ S7 give the severe constraint as seen in Figure 3. This CP violation phase 
also contributes on the CP-violating asymmetry of the non-leptonic decays of B® and B° s 
mesons. 

Let us discuss Sf, which is the measure of the CP-violating asymmetry, for B® — > 
J/ijjK s , 4>K S , rj'K . As discussed in Section 2, these <S/s are predicted to be same ones in 
the SM. On the other hand, if the squark flavor mixing contributes to the decay process at 
the one-loop level, these asymmetries are different from among as seen in Eq. ffl9|) . 

Since the phase 9^ contributes to A b (^ SJ , S^Kg and S v >k° of B® decays. We expect the 




Figure 3: A b ^ sl versus 6™ for \(8% R ) 23 \ = 
10 _3 (blue) and 10 _4 (orange), where the 
pink lines denote the experimental upper 
and lower bounds at la level. 

1.0 1 



0.8 



0.4 
0.2 




-0.06 -0.04 -0.02 0.00 0.02 0.04 



Figure 5: S^Kg of B® versus A^^ for 
|(^ fl ) 23 | = 10~ 3 (blue) and l(r 4 (orange). 



Figure 4: The allowed region on 1(5^)231- 
6*23^ plane. The experimental data at 90% 
C.L. of the branching ratio and Aqp^ 7 are 
taken account. 



1.0 




-0.06 -0.04 -0.02 0.00 0.02 0.04 



Figure 6: S^k of B® versus A 1 ^ 1 for 
\(5% R ) 23 \ = 10~ 3 (blue) and 10" 4 (orange). 



S 



correlations among them. We fix 1(5^)231 = 10 _4 (orange) and 10 _3 (blue) for typical values 
in the following calculations, we show the predicted regions on the A^^-S^Ks and A^^ 7 - 
S^k® planes in Figures 5 and 6, respectively. The experimental data is denoted by red lines 
at 90% C.L.. We also present the predicted region on the S^Ks'S^k plane in Figure 7, the 
slant dashed line denotes the SM prediction Sj/^k s = S<t>K s = SrfK, where the observed value 
Sj/^k s = 0.671 ± 0.023 is put. The reduction of the experimental error of A 1 ^^ will give 
us severe predictions for S^Ks an d <S v 'k°- It is noticed that this predicted region is different 
from the one in the previous work [12], where (5^)23 is neglected and /itan/3 = 0(10 TeV). 

As seen in Figures 5, 6 and 7, the reduction of the experimental errors will provide 
powerful tool to find the contribution of the squark flavor mixing. At last, we show the 
correlation between ^j^ 7 anc ^ °f the decay in Figure 8. We expect the observation 
of the CP violation in B° s <p<p at LHCb. 




0.0 0.2 0.4 0.6 0.8 1.0 -0.06 -0.04 -0.02 0.00 0.02 0.04 



•V s A C p Sy 

Figure 7: Predeicted region on the S^Kg- Figure 8: Predicted asymmetry in the 

S v , K o plane, where the slant dashed line B° s decay versus A b ^, sl . 
denotes the SM prediction. 



5 Summary 

The CP violation of the neutral B meson is the important phenomenon to search for new 
physics. We have discussed the contribution of the squark flavor mixing from (5^)23 and 
(5f L )23 on the direct CP violation of the b — > 57 decay and the CP- violating asymmetry in 
the non-leptonic decays of B® and B° s mesons. 

The magnitude of | (5^)23) is bounded by the branching ratio of b — > 57 with the con- 
straint of A^^ . The predicted branching ratio becomes inconsistent with the experimental 
data at \(5^ R ) 2 3\ > 5.5 x 10~ 3 . We have obtained the allowed region on the K^^sh^^f 
plane. While the | 23 1 is cut at 5.5 x 10 -3 , CP-violating phase Q^z is severely con- 
strained at 1(5^)23! — 2 x 10~ 3 . This CP-violating phase also contribute to the CP-violating 
asymmetry in the non-leptonic decays of B® and B° a mesons. 

We have predicted the correlation among A^~p S7 and Sf of the B^ and B° s decays. These 
CP-violating asymmetries could deviate from the SM predictions. 
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Since we suppose rather small /itan/3, 0(1 TeV), the contribution from LL(RR) compo- 
nents are minor in these CP- violating asymmetries. In this case, the new physics contribution 
is minor in M\ 2 of B q — B q (q = d, s) mixing since 1(5^)231 is at most O(10~ 3 ). This result 
is consistent with the recent result of the CP violations at LHCb as discussed in Section 2. 

In the near future, the precise data of the direct CP violation of b — > and CP- violating 
asymmetries in the non-leptonic decays of B® and B° s mesons give us the crucial test for our 
framework of the squark flavor mixing. 
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